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Abstract. We have studied systematically the influence of particle-hole 
symmetric and asymmetric kinetic terms on the ordered phases that we may 
observe competing or coexisting in a tetragonal system. We show that there 
are precise patterns of triplets of ordered phases that are accessible (i.e. it is 
impossible to observe two of them without the third one) . We found a systematic 
way to predict these patterns of states and tested it by identifying at least 16 
different patterns of three order parameters that necessarily coexist in the presence 
of the kinetic terms. We show that there are two types of general equations 
governing the competition of all these triplets of order parameters and we provide 
them. 
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1. Introduction 

Almost all important functional materials undergo a pleiad of phases that under 
certain conditions may coexist. Controlling the parameters leading to coexistence 
is of primary importance as it can provide access to new intriguing phenomena and 
functionalities. Already in the early eighties, systematic theoretical investigations of 
the coexistence of two ordered electronic phases appeared [1, 2, 3, 4, 5, 6, 7, 8, 9], 
motivated essentially by the general problem of antiferromagnetic superconductivity 
that emerged in organic superconductors and heavy fermions. Numerous theoretical 
studies of the coexistence and competition of two phases continue to appear as the 
number of such experimental paradigms multiplies. 

It has been a general conclusion from all the above studies that we must take into 
consideration on the same footing the two order parameters (OPs) that compete and 
may eventually coexist, otherwise we miss qualitatively new phenomena associated 
with this competition. However, as we shall show below, additional order parameters 
may coexist as well, and the need to include them is equally important. In fact, 
we will show that the particle-hole symmetric and asymmetric kinetic terms (KTs) 
of the hamiltonian, impose patterns of three order parameters (or triplets of order 
parameters) that are unavoidable. Whenever two of the order parameters coexist the 
third one appears as well. Therefore, we must necessarily consider all three order 
parameters simultaneously. 

We have considered a tetragonal tight binding system and we have studied 16 
cases of phase coexistence involving various types of ferromagnetism, density waves 
[10, 11, 12, 13, 14] and superconductivity [15]. In the case of a tetragonal lattice, the 
possible OPs that can be observed along with the kinetic terms are 63, generating 
an SU{8) Lie algebra [16]. The OPs we have studied, were chosen among those of 
an 5*0(8) subalgcbra [17, 18] that describes only even parity order parameters and 
is relevant for the nearly half-filled case. We have observed that the kinetic terms 
impose phases that do not initially exist in the Hamiltonian. These phases ought to 
have been already included in the Hamiltonian from the very beginning in order to 
study the system consistently. We have performed this procedure in all cases that 
we have studied and we obtained in all cases the self-consistence equations of all the 
order parameters involving the induced order. We observed that the two initial order 
parameters, the induced order and the mixing kinetic term satisfy a system of self- 
consistence equations that entangles their dynamics. They constitute closed sets of 
order parameters that need to be treated on the same footing. 

Through detailed examination of a number of systems with many phases 
emphasizing on the role of particle-hole symmetric and asymmetric kinetic terms, 
we managed to extract a simple empirical rule which helped us to predict patterns of 
OPs that coexist when the kinetic terms are properly taken into account. According to 
this rule the m,a,trix product of the matrix representations of the two initially coexisting 
order parameters and, the mixing kinetic term yields the matrix representation of the 
induced phase. Equivalently, the matrix product of the three involved order parameters 
yields the mixing KT matrix. In all cases that were selected according to this rule the 
predicted phase coexistence was confirmed. 

We have to remark that although the above rule is quite expected if we want a 
non zero mean value for the induced phase, there is no way to be certain that if this 
rule stands for a specific set of order parameters then we must obtain the coexistence 
described above unless we do the calculations. Indeed, only by considering on a 
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suitable spinor formalism the relevant triplets of order parameters within a BCS like 

mean field approach and extracting sclf-consistcnt gap equations wc were able to 
identify definitely that these states appear altogether as an unavoidable pattern. 

2. Results 

We consider all even parity OPs that are possible in a tetragonal system and may 
be relevant for discussing a number of heavy fermion materials as well as high-Tc 
cuprates. To describe in a unified way the coexistence of various OPs we need to 
introduce an eight component spinor formalism. We introduce the following spinor 

^fc = ( "^It "fci '^fc+QT '^fe+Qi "^-fcT '^-fci ^-k-Q^ Q^-fc-Qi ) (1) 

where afeg/ct^.^ are the destruction/ creation operators of an electron of momentum k in 
the Reduced Brillouin zone and spin projection s =t, i- This enlargement of the spinor 
space allows the simultaneous description of ferromagnetism, zone center (zero Cooper- 
pair momentum) and staggered (finite Cooper-pair momentum) superconductivity, 
charge and spin density waves. The density waves and the staggered superconductivity 
are characterized by the wave- vector Q = {it, it) which is the best nesting vector close 
to half-filling. To work in this eight dimensional spinor space we consider a base 
formed by the Kronecker products of the unit matrix and three of the usual Pauli 
matrices , pj, where i, j, fc = 1, 2, 3. 

The possible order parameters arising from the preceding spinor theory are 
4 X 4 X 4 = 64. If we demand that our Hamiltonian is traceless we are left with 
63 order parameters (including K.T.) that constitute the generators of an SU{8) 
spectrum generating algebra [16]. In the case of tetragonal systems close to half- 
filling, equivalence of the Brillouin zone points (7r,0) = ((—• tt, 0)) and (0,7r) = (0, —w) 
imposes that the order parameters have even parity. The OPs satisfying this constraint 
are 28 (including K.T.) and form an SO{8) spectrum generating algebra [17, 18]. The 
OPs that we have considered in this study were chosen among those 28 identified in 
Table 1 with their symbols adopted here. 

Table 1. The 28 OPs that form an SO{8) spectrum generating algebra and would 
be accessible in a tetragonal system close to half filling. In the next sections we 
demonstrate that particle-hole asymmetric and symmetric kinetic terms impose 
various patterns of triplets of the following OPs. 



Or-der Fammcler Type 



7 


nearest neighbours hopping term 


5 


next nearest neighbours hopping term 


^ x,y,z 


ferromagnet along x,y,z-axis 


A 


d-wave ferromagnet along x,y,z-axis 


w 


charge density wave 


Jc 


orbital anti-ferromagnet 




spin density wave along x,y,z-axis 


x,y,z 


spin nematic along x,y,z-axis 


A. 


s-wave SC (q = 0) 




d-wave SC (q = 0) 


V 


s-wave SC (q = Q) 




d-wave SC along x,y,z-axis (q = Q) 
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We note that in Table 1 there are 16 OPs corresponding to particle-hole 

condensates including the KTs and 12 superconducting states (including their complex 
conjugates). Moreover, 8 of the 12 superconducting OPs represent staggered 
SC in which the pairs have a finite total center-of-mass momentum {q = Q) 
bearing similarities to the Fulde-Ferrel states [19]. These quite exotic states are 
superconducting states with modulated superfluid density and as we will show below, 
they should play a crucial role in any antiferromagnetic SC state. 

We report here 16 different patterns of OPs that are imposed by the particle- 
hole symmetric and asymmetric kinetic terms. They can be classified into 3 different 
types of OP mixing that according to their properties can be merged into two general 
groups. In all these cases we present the typical system of self-consistence equations 
that provide the OPs and we identify the kinetic terms that are responsible for the 
OPs mixing. 

2.1. First type of OPs mixing 

In this first case, we consider that the Hamiltonian consists of the two kinetic terms 
and three order parameters. These order parameters have been chosen according to 
the empirical ride mentioned in the introduction i.e. the matrix product of the three 
order parameters yields the kinetic term that causes their mixing. In Table 2 we 
present the different combinations that fall into this class. 

Table 2. Triplets of order parameters that form patterns imposed by the kinetic 
terms mentioned in the last column. In all cases we have the same system of 
self-consistence equations for the OPs, provided we replace the corresponding 
OPs of the same column. 



OP 1 


OP 2 


OP 3 


mixing KT 


M, 






5 


W 


A. 


V 


5 




A. 




5 


n. 


A. 




7 


V 


Ad 


W 


7 








7 



We consider explicitly the first combination of the preceding table in order to 
demonstrate the general equations governing the phase coexistence and competition 
in the above patterns. In fact, the self-consistence equations that result are the same 
for all six patterns provided we replace the corresponding OPs that are in the same 
column. The Hamiltonian corresponding to the first case is given by the relation 

fe 

where we have suppressed the momentum index k and the Kronecker product's symbol 
(g). The energy eigenvalues are 



E± = Jm| + 72 + ^2 + n2 + A2 ± 2^ (M| + 72)52 - 25M,n^A^ + {A^ + 72)02 (3) 
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The self-consistence equations of the order parameters are the following 



M(El-El-45^)+46AJl^ ^ f E^\] 

eT^ — [w)] 

A^(El-El-ml)+4dM^IL^ ^ /E^\\ 

e7^ ^ j I (5) 

Il(El-El-4j^-4Al)+4SA^M^ ^ / E\\ 

^ tanh j I (6) 

where the OPs and KTs depend on k' . The mixing role of the kinetic term is explicit 
already in the form of the equation. In fact, the usual BCS equations for each one of 
the order parameters are expected to have the general form 



fe' 

where / is a function of the energy dispersion E and temperature T. Each BCS 
equation supports solutions of zero and non-zero order parameter, depending on the 
temperature. On the other hand, in our case each OP self-consistence equation has 
the general form 

\ = E^fe,fe' {fiE,„T)A^,+g{E,„T)m^,A^,B^,} (8) 
fe' 

where /, g are function of the energy dispersion and temperature, rUf^ is the mixing 

kinetic term and A^. and Bj^ are the other two OPs. We observe that a solution of 
zero order parameter is not possible unless the mixing term or one at least of the other 
order parameters is also zero. This suggests that in the presence of the mixing 
kinetic term we cannot have two order parameters without the third. The 
three order parameters and the mixing term constitute a group that must be treated as 
an independent subsystem on the same footing. 

It is interesting to obtain the sclf-consistcncc equations when the kinetic term 
that does not contribute to the mixing is set to zero. In this case, the eigenenergies 
obtain the form 



s+ = V(^.+<5)2 + (A,-n,)2 

E_ = ^{M,-Sr + {A, + U,)^ 



(9) 

(10) 
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tanh 



2T 

, fE- 
tann —— 
\2T 



(11) 



^(M,+5)2 + (A,-n,)2 
/(M, -5)2 + (A^ + n,)2 

2.2. Second type of OPs mixing 

The second case involves a different coexistence pattern involving once again three 
order parameters and a mixing kinetic term. The following table contains the 
combinations belonging in this class: 

Table 3. Same as in Table 2 



OP 1 


OP 2 


OP 3 


mixing KT 






Tly 


5 


A. 


ri 




5 








7 


A, 


w 




7 



Once again we present the typical results of one of these cases, specifically we consider 
the first. The Hamiltonian is 

^ = 51 (^^3 P3 + ^^3 - ^y-^l Pi - ^x-^l P2<^3 ' ^3 ^^3) *fe (12) 
fc 

The corresponding quasi-particle poles are 



Finally, we obtain the self-consistence equations 



(13) 
(14) 




F, -5 



+ 



(f, -5) + (n^ + Hj,) 



tanh 



tanh 



2T 



E4 



2T 



tanh 



— tanh 



E-+ 
2T 



2T 



(15) 
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+ 







V V 


+ IlyY +(F^ 























tanh 



tanh 



tanh 



tanh 



E. 



++ 



2T 



2T 



E- 



++ 



2T 



E. 



2T 



tanh 



tanh 



tanh 



tanh 



E^ 



2T 



E+- 
2T 



ki6) 



2T 



E. 



2T 



ki7) 



As one can observe, the above equations do not have the general form (8) in which the 
mixing role of the relevant kinetic term is explicit as in the previous case. However, by 
performing a Taylor expansion with respect to the order parameters up to quadratic 
order terms we can show that such a relation does exist. These self-consistence 
equations imply once again that if the mixing term is present we cannot have the 
two order parameters without the third. 

Indeed, let us consider (15) supposing that is absent from the initial 
Hamiltonian. Then we set F^_ = in the right side of (15) and this equation will 
now provide the induced part of F^ : 



I induced 



E 



+- 



(52 + (n^ + n^y 
s 



tanh f tanh 



E' 

tanh ( — 
2T 



■ tanh 



E'_ 



2T 



^j5^+{ll,-Uyy 

where we have introduced the new energy dispersions E', by setting F^ — 0: 



(18) 



E'^_=^±^S^+(ll,-Ily) 



(19) 
(20) 



A zero induced F^ term is expected if one of the two following conditions holds. 
On one hand, we may have 5 = 0, in which case we confirm that we will not have 
induced F^ if the mixing kinetic term vanishes. On the other hand, we may have 
E±+ = and E±- = -Bq:-. This last condition can be realized only when 11^ = 
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or Ily = 0. Consequently we conclude that if the mixing term is present and two order 
parameters are non zero we have an induced order . i.e. the three phases coexist. 

It is instructive to derive here as well, the self-consistence equations when the 
irrelevant, to the mixing, kinetic term vanishes. 



E±+ = El-^ = ±^ (^F, - 5) + (n, + Hy) (21) 

E±- = Elz° = ±^(f. +<5)'+ (n,-n^)' (22) 

We observe that E+± = —E-^. This equality simphfies the self-consistence equations. 
For example we have 




Quite remarkably, we have encountered the same form of self-consistence equation in 
Section 2.1 when the non mixing kinetic term was set to zero. This common feature 
reveals that these two cases share the same mixing "mechanism", constituting specific 
examples of a more general coexistence pattern. 

The next pattern that we shall discuss here is the coexistence of two specific 
phases, s-wavc and d-wavc SC OPs in the presence of the two kinetic terms, where the 
kinetic terms play both the role of the mixing terms and the OPs at the same time. 
As far as the form of the equations that we derive, they belong to the same general 
coexistence pattern like the one reported just above. The Hamiltonian is 

^ = 13 *fe (^^3 P3 + '5t3 - A^r2 P3a2 - A,T2 *fe (24) 
k 

The poles of the Green's function are 



£± = '^(a,±A,) +{^±5f (25) 

It is evident that they have the form of the previous cases Section 2.1 and Section 2.2 
whc'ii wc set the irrelevant kinetic terms equal to zero. The self-consistence equations 
obey the same rule 



^1-^l^^kk'\ r- Z^^''^\2T 
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fe' 



+ 



A, -A, 



J(a 


— A 


J +{'J 0} 




A,^ 


^A, 




+ A, 






A, 


-A, 




-A, 





tanh ( 



tanh 



2T 



tanh ( ^ 



(26) 



2.3. Third type of OPs mixing 

The next case we consider has distinct properties from the preceding encountered in 
Section 2.1 and Section 2.2. We have found that the following combinations all have 
the same coexistence pattern. 



Table 4. Distinct type of mixing compared to the one related to Section 2.1 and 
Section 2.2. The following combinations obey the same system of self-consistence 

equations. 



OP 1 


OP 2 


OP 3 


mixing KT 




W 




5 




Jc 


Ji 


5 




•^y 




S 


F, 


V 




7 



The example in this type of mixing is the first of Table 4, which is described by the 
Hamiltonian 

^ = *I (^^3P3 + ^7-3 - M^T3PiC73 - WtsPi - F^T^G^ (27) 
k 

with corresponding eigenenergies 



e±_ = {f,t5)-J(m,±w) +72 



(28) 
(29) 



We observe that the structure of the poles are different from the ones found in Sections 
2.1 and 2.2. The self- consistence equations are given from the following relations 



P^=\Y. {tanh (^) + tanh {^^) + tanh (^) + tanh (^) } 

(30) 
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We have to remark that Equations (31) and (32) have similarities with the results of 

the previous sections. Though, Equation (30) is totally different. Close observation of 
the equation and the eigenenergies, shows that great simplification occurs when 7 = 0. 
In this case we have 



E^^ = {F,tS) + {M,±W) (33) 
E^_ = {F,tS)-{M,±W) (34) 
and as far as the self-consistence equations are concerned 




(37) 

According to what he have been taught from the previous sections, the symmetry 
the above equations present, implies that we have reached to a triplet of order 
parameters that necessarily coexist in the presence of the corresponding mixing kinetic 
term. This can be shown as follows. Any of these order parameters can be zero only 
if the spectrum is particle-hole symmetric. This occurs only when two out of the four 
terms are zero. Consequently if three of these terms are non zero the fourth will be 
non zero, too. 
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The final case we present, has two order parameters and the two kinetic terms. 
As wc shall see it belongs to the same coexistence pattern of the above cases. The 
Hamiltonian is 

^ = E *I {^^sPs + + ^.7-3^3^1 + F^T^a^ (38) 

k 

The eigenenergies are equal to 

£;±+ = (A,.±Fj + (<5±7) 
i?±_ = (yl,±Fj-(<5±7) 
The corresponding self-consistence equations are 




(42) 

3. Discussion 

Having studied 16 cases of coexisting OPs in the presence of the particlc-holc 
symmetric and asymmetric KTs, we have found 3 different types of coexistence 
patterns. The 11 cases studied in Sections 2.1 and 2.2 seem to have the same type 
of coexistence "mechanism" belonging to a more general Coexistence Schem,e. On 
the other hand the 5 cases that we presented in Section 2.3 originate by a distinct 
coexistence "mechanism". Consequently, we conclude that the 16 cases we 
studied merge into two general Coexistence Schemes as in Table 5. 

Moreover, we have found that special coexistence patterns can be observed even 
when he have only 2 OPs. In this case the kinetic terms play a dual role. They 
behave as the mixing KTs and the members of triplets of coexisting OPs. In that 
case, with both kinetic terms present we cannot observe one of the two OPs without 
the second one. For example, in the presence of both kinetic terms, d-wave SC 
coexists with s-wave SC ! This shows that the generation of these patterns is not a 
special property owned by the KTs. Consequently, these Coexistence Schemes 
originate due to more general relations that are satisfied by quartets of 
the SU{S) generators. We expect that other terms, apart from the KTs could 
play the role of the mixing terms and produce different type of coexistence patterns. 
We will present elsewhere a complete account of all patterns of coexisting states that 
correspond to the above mentioned quartets. 

Finally, apart from the general conclusions that we inferred about phase 
coexistence, we obtained valuable results concerning specific coexisting triplets of 
OPs that may correspond to the physical situation in numerous correlated systems of 
interest. Particularly, we have observed that: 



(39) 
(40) 
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Table 5. The two distinct Coexistence Schemes 



Scheme 


OP 1 


OP 2 


OP 3 


mixing KT 


I 


M 




n 


s 


I 


W 


A, 




s 


I 


J" 






s 


I 




A, 




'Y 


I 




^d 


w 


7 


I 


Uy 


^d 


,/= 

y 


7 


I 


Fz 




Hy 


S 


I 


Az 


v 


n. 


s 


I 


Fy 


J'. 


Mz 


7 


I 


Az 


w 


Mz 


7 


I 


A. 




5 


7 


II 


Fz 


w 


Mz 


(5 


II 


Fz 


Jc 


Ji 


S 


II 


Ax 


Jy 




5 


II 


Fz 






7 


II 


J:\, 


X 


<) 





• Density waves, zone-center superconductivity {q = o) and staggered super- 
conductivity (q = Q) constitute a triplet of OPs that necessarily coexist in 
the presence of the KTs. Such an observation is of general relevance for all 
antiferromagnetic superconductors, a category of materials that includes organics, 
heavy fermions, high- Tc cuprates etc. 

• Ferromagnetism, charge density waves and spin density waves constitute another 
triplet of OPs that necessarily coexist in the presence of the asymmetric KT. This 
observation has already been reported before and shown to be related with the 

colossal magnctoresistancc phenomenon [20]. 

• s-wave and d-wave superconducting OPs always coexist in the combined presence 
of the symmetric and asymmetric KTs. Needless to note that high- Tc cuprates 

as well as numerous heavy fermion systems arc believed to be d-wave SC. Our 
observations imply that a pure d-wave SC state is an oversimplification. 

• s-wave and d-wave ferromagnetic OPs always coexist in the combined presence of 
the symmetric and asymmetric KTs. The implications of this observation need 
to be investigated. 
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